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ABSTRACT 


Expressions  for  the  vibratory  pressure  field  produced  by  an 
operating  counterrotating  propeller  system  are  developed  In  terms  of 
first  and  second  blade  harmonics  of  the  individual  propellers.  The 
coefficients  of  the  first  and  second  blade  harmonics  of  a  3-bladed 
propeller  are  expressed  in  closed  form  in  terms  of  complete  elliptic 
integrals..  It  is  found  that  the  pressure  signal  due  to  a  counter¬ 
rotating  propeller  can  be  obtained  from  the  summation  of  the  pressure 
fields  of  its  two  components  as  though  the  two  propellers  operate  as 
separate  units,  since  their  mutual  interference  has  been  found  to 
contribute  little  to  the  total  vibratory  pressure  signal  in  the  near 
field.  This  study  indicates  that  the  counterrotating  propeller  system 
has  vibratory  characteristics  much  superior  to  an  equivalent  single 
propeller  having  the  number  of  blades  of  one  of  its  components,  delivering 
the  same  thrust  at  the  same  RPM  at  optimum  conditions.  A  study  of  the 
interference  effects  by  means  of  a  two-dimensional  analogue  indicates 
that  the  effect  of  the  forward  propeller  on  the  aft  is  much  greater 
than  the  effect  of  the  aft  propeller  on  the  forward  propeller. 


INTRODUCTION 


Counterrotating  propellers  are  an  important  means  of  implementing 
high-powered  propulsive  units.  It  has  been  generally  accepted  that,  with 
this  propeller  system,  propulsion  losses  are  minimum  and  therefore 
efficiency  is  high.  However,  because  of  mechanical  complications  this 
type  of  propulsion  has  found  little  practical  application*  it  has  been 
adopted  chiefly  for  torpedoes.  It  is  the  purpose  of  this  investigation 
to  consider  the  vibration-producing  characteristics  of  counterrotating 
propellers  and  to  compare  them  with  those  of  a  single  propeller  with  the 
expectation  that  counterrotating  propellers  may  show  advantages  in 
reducing  vibration. 

The  pressure  field  generated  by  an  operating  propeller  is  mainly 
responsible  for  the  vibrations  experienced  by  neighboring  boundaries. 

In  view  of  present  technical  developments  with  continuously  increasing 
propulsive  forces,  the  vibratory  motions  due  to  the  pulsating  pressure 
field  may  well  be  expected  to  Increase  to  a  degree  which  will  require 
either  a  drastic  change  in  the  stern  configuration  or  a  different 
propulsion  device.  Although  other  considerations  such  as  avoidance  of 
cavitation  and  diminution  of  noise  will  undoubtedly  weigh  more  heavily 
in  the  selection  of  a  new  design  than  the  vibration — causing  qualities, 
nevertheless  it  would  appear  important  in  future  design  work  to  evaluate 
these  characteristics  for  the  two  configurations  which  are  feasible 
changes  from  the  usual  propeller  arrangement,  viz,  counterrotating 
propellers  and  shrouded  or  "nozzle"  propellers. 

Breslln  (l)  has  made  a  theoretical  study  of  the  pressure  field 
near  an  operating  ship  propeller  in  the  open  water  condition.  By  using 
vortex  line  theory  and  representing  the  propeller  blade  by  three  types 
of  vorticity,  a  hub  vortex  streaming  aft  along  the  x-axls,  a  bound 
vortex  fixed  at  a  position  of  the  blade  axis  and  a  tip  vortex  streaming 
aft  in  the  form  of  a  helix,  he  determined  the  instantaneous  pressure 
at  any  point  in  the  field. 

It  was  found  that  the  total  pressure  at  a  point  in  the  field  is 


purely  convective  pressure  ('^3^/9x)  since  the  impulsive  part  of  pressure 
due  to  the  time  rate  of  change  of  blade  position  is  annulled  by  part  of 
the  convective  pressure.  It  may  be  stated,  therefore,  that  the  pressure 
can  be  investigated  as  a  steady-state  phenomenon  depending  on  the  blade 
position.  This  fact  makes  for  considerable  simplification  in  the  present 
investigation.  In  References  2  and  3,  the  blade  frequency  pressure  and 
velocity  fields  around  a  marine  propeller  due  to  loading  and  thickness 
effects  have  been  determined  in  closed  form  in  terms  of  complete 
elliptic  integrals.  If  the  mutual  interference  of  the  two  propellers  is 
neglected  and  the  results  of  References  2  and  3  are  used,  the  pressure 
field  around  two  counterrotating  propellers  will  be  determined  as  a 
function  of  blade  position.  In  contrast  to  the  case  of  a  single  propeller 
where  only  the  m-blade  frequency  need  be  considered,  in  the  case  of  two 
counterrotating  m-bladed  propellers  both  the  m  and  the  2m  blade 
frequencies  and  — )  are  of  practical  significance. 

2lT  TT 

In  the  first  part  of  the  present  investigation  expressions  for 
the  fluctuating  pressure  signal  due  :o  the  counterrotating  propellers 
have  been  developed  as  though  the  two  propellers  operate  as  separate 
units  without  any  interference  between  them.  Expressions  for  the  first 
and  second  blade  harmonics  of  the  pressure  signal  have  been  evolved 
for  the  case  of  a  pair  of  3-bladed  propellers  in  terms  of  complete 
elliptic  integrals. 

The  second  part  of  this  investigation  is  concerned .with  the  mutual 
interference  of  blades  as  they  approach  and  pass  each  other.  The  inter¬ 
ference  leads  to  a  time-dependent  increment  of  lift  (thrust)  which 
averaged  over  time  remains  close  to  zero  but  instantaneously  may  be 
large  although  of  short  duration.  This  interference  effect  is  studied 
under  two-dimensional  flow  conditions  which  are  more  amenable  to 
mathematical  analysis  than  the  three-dimensional  flow.  The  magnitude 
of  the  interference  effect  on  the  pressure  in  an  actual  three-dimensional 
case  will  be  much  smaller. 

The  third  part  deals  with  the  evaluation  of  the  pressure  field 
in  the  presence  of  3-0-3  counterrotating  propellers.  (The  first  digit 
in  the  designation  denotes  the  number  of  blades  in  the  forward  component, 
the  second  the  number  of  countervanes  and  the  third  the  number  of  blades 


in  the  after  component.)  The  expressions  evolved  in  the  analysis  are 
used  to  determine  the  coefficients  of  the  3rd  and  6th  harmonics  of  the 
pressure  signal  due  to  each  3-bladed  propeller  separately  at  6  different 
stations,  up  to  an  axial  distance  of  0.8  of  the  propeller  radius,  for  four 
different  tip  clearances.  By  combining  the  corresponding  Fourier  coefficients 
the  pressure  signal  in  the  field  of  the  3-0-3  propeller  system  is  determined. 
In  addition,  the  case  of  a  single  3-bladed  propeller  producing  equivalent 
thrust  and  having  the  same  RPM  as  the  3-0-3  system  is  treated  and  its 
pressure  field  compared  with  that  of  the  counterrotating  propellers. 

This  study  is  prepared  under  Contract  Nonr  263(16),  of  the  Bureau  of  Ships 
Fundamental  Hydrodynamics  Research  Program,  Project  S-R009-01-01,  and 
technically  administered  by  the  David  Taylor  Model  Basin. 


ANALYSIS 


Part  I 

Simplified  Theory  for  the  Oscillatory  Pressure  Field 
About  Counterrotating  Propellers 

In  this  propeller  system  two  coaxial  propellers  are  situated  a 
short  distance  apart  and  are  rotated  in  opposite  directions.  Due  to  the 
short  axial  separation  there  will  be  mutual  Interference  between  the 
two  propellers.  In  this  part  of  the  Investigation,  however,  the  blade 
Interference  effects  will  be  neglected  and  the  pressure  field  will  be 
evaluated  as  If  the  two  propellers  Influenced  the  field  separately. 

Breslln  (1)  has  shown  by  vortex  line  theory  and  also  by  use  of 
doublet  distributions  that  the  linearized  pressure  charge  In  an  other- 
wise  unbounded  fluid  due  to  an  operating  m-bladed  propeller  Is  the  sum 
of  thrust-produced  pressure  and  torque-produced  pressure  as  followst 
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T'  ■  thrust  loading  coefficient  •  T/vd  where  T  Is  thrust  and 
b  propeller  radius 
m  ■  number  of  blades 
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J  ■  — fr  ■  advance  ratio,  where  v  Is  the  velocity  of  advance 
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n  -  r.p.s.  and  D  propeller  diameter 

/  2  2  2  ^  sl/2 

R  ■  (x  +r  •♦•s  -  2rs  cosa; 

0 

x,r,Y  “  cylindrical  coordinates  of  the  field  point  In  the  fluid  at 
which  the  pressure  signal  Is  calculated* 

■  0  _  Y  ,  where  9  is  the  blade  position  angle. 

s  Is  dimensionless  blade  span. 

The  foregoing  expressions  (l)  and  (2)  are  referred  to  a  coordinate 
system  with  the  origin  at  the  propeller  center  (Fig.  l)  and  all  the  linear 
dimensions  are  modtiples  of  the  propeller  radius  b  . 

In  the  counterrotating  propeller  configuration,  the  origin  of 
the  coordinate  system  is  taken  on  the  axis  of  rotation  midway  between 
the  propeller  mid-planes,  i.e.  at  a  distance  ^  from  the  forward 
and  aft  propeller,  respectively.  The  x-axis  is  taken  along  the  axis 
of  rotation  opposite  In  direction  to  the  main  stream  flow  (Fig,  2), 

Since  Eqs,  1  and  2  have  been  developed  by  making  use  of  the 
results  of  the  linearized  theory,  the  instantaneous  pressure  signal 
of  the  counterrotating  propeller  configuration  will  be  given  by 

P  -  P^2))  ^  (3) 


where  the  forward  and  the  after  propeller  are  designated  by  superscript 

12  12 

1  and  2,  respectively,  and  P^*  and  Pp’  are  as  given  in  Eqs.  (1)  and  (2) 

with  the  only  difference  that  x  should  be  replaced  by  x  -  x^  for 
the  forward  propeller  and  by  x  +  x^  for  the  after  propeller. 

An  important  aspect  of  the  pressure  field  is  revealed  by  seeking 
the  contribution  at  m  and  2m  frequencies.  The  total  pressure  signal 
at  any  given  point  of  the  field  will  be  expressed  in  terms  of  Fourier 
coefficients  as 
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where  0^-Y  ■  the  Instantaneous  angular  position  of  the  blade  relative 

to  the  station  where  the  pressure  signal  is  determined.  In  the  first 
summation  over  1,1  and  2  correspond  to  the  forward  and  aft  propeller. 
The  second  summation  over  k  gives  the  contribution  of  all  harmonics 
and  the  third  summation  from  n  ■  1  to  m  -  1  arises  through  the 
Inclusion  of  the  contributions  of  the  other  m  -  1  blades  which  are 
spaced  at  angular  Intervals  of  ~  ~ 

It  Is  shown  in  Appendix  A  that 
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for  k  not  an  Integer  multiple 
of  m  . 
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for  k  an  Integer  multiple 
of  m  . 


For  the  case  of  k  a  non- Integer  multiple  of  m  ,  P  ■  0  so 
that  there  is  no  contribution  to  the  pressure  signal  from  the  harmonics 
which  are  not  integer  multiples  of  the  blade  frequencies.  The  only 
terms  which  contrlbyte  to  the  pressure  signal  are  harmonics  of  m  Oj 
or  higher  multiples  of  m  a  »  which  are  known  as  the  harmonics  of 


the  blade  frequencies.  Therefore,  the  total  pressure  signal  will  be 
given  as 
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When  only  the  first  two  blade  frequencies  are  considered. 
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The  angular  blade  positions  and  for  the  two  components  of 

the  counterrotating  configuration  change  with  time  according  to  the 
relationship 


since  both  are  rotating  with  the  same  angular  velocity  u  around  the  same 
axis  but  In  opposite  directions.  Then  the  angular  blade  position  of  the 
after  propeller  Is  given  In  terms  of  the  forward  by 


2(»)t>c 


where  c  Is  a  constant  depend5.ng  on  the  Initial  conditions.  Under  the 

assumption  that  at  t  ■  0  ,  0  and  —  ,  the  above  relation 

m 

becomes 
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It  Is  of  Interest  to  examine  the  situation  at  two  Instants  of 
time^  when  the  blades  are  In  alignment  (in  phase)  and  v^en  they  are 
farthest  from  alignment  (out-of-phase).  The  In^phase  condition  will 


exist  when 
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and  the  out-of-phase  condition  will  be  given  when 
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i.e,  when  the  blades  of  the  after  propeller  bisect  the  angles  between 
the  blades  of  the  forward  propeller.  For  the  in-phase  condition,  the 
total  pressure  signal  will  be  given  (see  Appendix  B)  by 
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of  the 

Fourier  series  representing  the  expressions  for  the  pressure  components 
of  Eq, (2),  For  the  out-of-phase  case,  the  total  pressure  signal  will 
be  given  (see  Appendix  B)  by 
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Therefore  for  all  instants  at  which  the  blades  are  in  coincidence  the 
pressure  signal  contains  m-blade  frequencies  as  well  as  2m-blade 
frequencies,  whereas  out-of-phase  it  seems  that  the  m-blade  signals 
tend  to  cancel  each  other  leaving  mainly  the  2m-frequencles .  As 
will  be  seen  later  on  from  the  calculations,  the  cancellation  will 
occur  at  large  distance  from  the  propeller  plane. 

In  order  to  calculate  the  pressure  signal  at  any  point  in  the 
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field  by  means  of  eq.  (9)  and  (lO)  it  is  necessary  to  evaluate  the 
Fourier  coefficients  for  the  first  and  second  blade  frequencies  in 
advance.  This  can  be  done  by  restricting  our  attention  to  the 
pressure  field  near  a  single  m-bladed  propeller.  Then  the  general 
expressions  which  have  been  developed  may  easily  be  transformed  to 
equivalent  expressions  for  the  description  of  the  pressure  field  due 
to  a  pair  of  counterrotating  propellers  at  a  distance  1.  from 
the  new  origin  (Fig.  2). 

As  demonstrated  in  Reference  2,  for  the  case  of  3-bladed 
propellers  the  thrust  contribution  contains  only  cosine  terms  and  the 
torque  contribution  only  sine  terms.  It  was  found  in  that  reference 
that  the  Fourier  coefficients  of  the  first  blade  harmonic  of  the 
pressure  signal  due  to  the  loading  are  given  by 
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for  the  cosine  term  and 
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for  the  sine  coefficient,  where 
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n(a^»k'),  r](p^  ',k')  are  complete  elliptic  integrals  of  the 
third  kind 

F(k*)  and  E(k')  are  complete  elliptic  integrals  of  the  first 
and  second  kind  respectively 


and  other  symbols  are  as  defined  in  eq,  (2), 

In  Appendix  C,  the  general  procedure  for  determining  the 
coefficients  of  the  second  blade  frequencies  of  the  pressure  signal 
due  to  a  3-bladed  propeller  and  the  evaluation  of  some  representative 
integrals  are  demonstrated.  The  results  aret 
for  the  cosine  coefficient. 
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and  for  the  corresponding  sine  coefficient. 
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Expressions  (ll)  (12)  (I4)and  (15 )have  been  used  together  with  eqs. 
(9)  and  (10)  to  evaluate  the  pressure  signal  of  the  3-0-3  counter¬ 
rotating  propeller  configuration.  At  x  »  0  ,  these  expressions  show 
signs  of  Indeterminacy  but  after  straight-forward  manipulation  In 
conjunction  with  a  limiting  procedure  they  are  reduced  toi 
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In  spite  of  the  fact  that  closed-form  expressions  have  been 
developed  for  the  fluctuating  pressure  signal  near  the  counterrotating 
propeller,  the  numerical  method  described  in  Reference  4  is  used  in  the 


present  work  since  a  coding  program  fur  the  evaluation  of  the  total 
Instantaneous  pressure  signal  was  available  for  the  ELECOM  100  computer. 

Hence,  as  in  Reference  4,  the  total  oscillatory  pressure  "signature" 
for  a  single  blade  is  plotted  for  the  angular  interval  -n  to  it  .  Since 
the  pressure  equation  is  linear  the  contribution  of  individual  blades 
can  be  superimposed  and  the  result  for  three  blades,  for  example,  can 
be  obtained  from  that  for  a  single  blade  by  summing  the  contributions 
of  the  one  blade  at  three  120°  Intervals.  The  resulting  composite 
curve  which  has  a  cycle  of  120°  for  three  blades,  can  be  analyzed 
harmonically  by  a  24-ordinate  numerical  procedure. 


Approximate  Evaluation  of  Pressure  Fluctuation 
Arising  from  the  Blade  Interference  Effects 

In  the  previous  section  the  vibratory  pressure  field  arising  from 
the  m-O-m  counterrotating  propeller  system  has  been  evaluated  as  if  the 
two  propellers  operate  separately.  It  is  known,  however,  that  because 
of  the  small  axial  distance  between  the  two  propellers  their  mutual 
interference  will  develop  a  pressure  field  which  will  give  rise  to 
a  time-dependent  lift.  It  is  believed  that  as  two  blade  sections 
approach  and  pass  each  other  in  relative  motion  this  time-dependent 
lift  remains  on  the  average  very  close  to  zero,  but  instantaneously 
it  may  be  large  although  of  short  duration.  The  interaction  phenomenon 
can  possibly  produce  strong  fluctuations  in  the  close  proximity  of  the 
propellers  and  therefore  it  should  be  included  in  the  analysis  of  the 
oscillatory  pressure  field  of  counterrotating  propellers.  A  study  of 
this  interference  phenomenon  as  a  three-dimensional  unsteady  flow 
problem  is  mathematically  untractab'e.  It  has  been  decided  therefore 
to  study  this  phenomenon  as  a  two-dimensional  problem  with  the  thought 
that  if  such  contributions  to  the  pressure  can  be  shown  to  be  small 
then  they  can  be  safely  assumed  as  small  in  frie  actual  three-dimensional 
case.  It  is  known  that  in  the  three-dimensional  mathematical  model 
(with  finite  span),  which  is  described  by  a  bound  vortex  line  and  free 
vortex  shed  from  the  tips  of  the  wing,  the  induced  velocity  field  is 
of  smaller  amplitude  than  is  that  of  the  corresponding  two-dimensional 
model  with  infinite  span  represented  by  a  lifting  line.  The  trailing 
vortex  reduces  the  angle  of  attack  of  the  incoming  flow,  and  consequently 
the  lift  and  the  generated  pressure  will  be  diminished. 

The  two-dimensional  problem  that  is  perhaps  most  closely 
related  to  the  one  posed  by  the  motion  of  the  propeller  blades  is 
that  presented  by  two  arrays  of  finite-chord  foil  sections  moving 
with  the  same  angular  velocity  Q  but  in  opposite  directions  as  shown 
in  Fig.  2b,  It  is,  however,  more  convenient  to  study  these  effects 
as  of  the  relative  motion  of  one  cascade  with  respect  to  the  other. 
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If  the  string  of  blades  of  the  forward  propeller  is  assumed  to  be 
fixed  in  space  the  cascade  of  the  after  propeller  blades  will  then 
be  moving  with  angular  velocity  2S2  .  This  sort  of  blade  arrange¬ 
ment  is  identical  with  that  encountered  in  the  stator  and  rotor  cascades 
of  axial-flow  machinery. 


In  order  to  preserve  as  many  as  possible  of  the  geometric  and 
kinematic  characteristics  of  the  three-dimensional  problem  in  this  two- 
dimensional  model,  so  that  the  latter  will  become  a  reasonably  useful 
analogue,  it  will  be  assumed  that  the  geometric  blade  characteristics 
are  those  existing  at  the  effective  radius  b^=  0,7b  of  the  propeller. 
The  blade  spacing  d  will  be  determined  as 


2TT(b  ) 

e  s ,r 


s,r 


m 


s.r 


where  m  is  the  number  of  blades  in  the  propeller  and  subscripts  s 
and  r  refer  to  the  stator  and.  rotor  respectively.  The  incoming  flow 
to  this  cascade  arrangement  has  a  velocity  V  in  the  axial  direction 
equivalent  to  the  speed  of  advance  and  a  tangential  velocity  equal  to 
the  product  of  effective  radius  and  angular  velocity. 

This  two-dimensional  analogue  has  been  studied  in  Reference  5 
as  a  useful  model  for  evaluation  of  the  induced  effects  on  the  blades 
of  a  typical  turbomachine  moving  through  a  nonuniform  field  disturbed 
by  the  stationary  blades.  The  results  of  the  theory  of  a  single  thin 
airfoil  in  non-uniform  motion  have  been  utilized  in  order  to  calculate 
the  resulting  unsteady  effects  on  an  infinite  cascade  of  blades  moving 
with  respect  to  a  stationary  row  of  blades  (stator).  Expressions  have 
been  developed  for  the  induced  velocities  and  the  lift  on  a  single 
rotor  blade  due  to  the  steady-state  circulation  distribution  of  each 
and  every  stator  blade  as  well  as  to  the  effect  of  the  vortlcltles 
shed  by  the  stator  blades  and  for  the  lift  on  a  single  stator  blade  due 
to  the  steady-state  circulation  distribution  of  each  rotor  blade.  Both 
flat  plate  and  elliptic  distr ibuticns  of  the  steady  circulation  are 
assumed  in  Reference  5, 


One  must  keep  in  mind,  however,  that  in  reality  the  stator- 


rotor  interference  problem  is  much  more  complicated  and  many  more 
factors  are  involved  than  are  treated  in  the  above  reference.  Each 
foil  is  influenced  by  l)  its  own  wake#  2)  the  variable  bound  vortices 
of  the  fellow  members  of  its  own  blade  row#  3)  their  wakes#  4)  the 
variable  bound  vortices  of  members  of  the  other  blade  row#  5)  their 
wakes  and  finally  6)  by  the  displacement  effects  of  the  blades#  since 
each  blade  is  of  finite  chord  and  thickness.  The  problem  can  be 
rigorously  treated  by  solving  a  pair  of  Integral  equations  which  state 
the  requirement  of  the  boundary  conditions  on  the  stationary  and  moving 
blades#  viz,  the  tangency  of  the  flow  along  the  moving  blades  and  zero 
velocity  across  the  stationary  blades.  The  solution  of  this  system  of 
Integral  equations  will  determine  the  circulation  distribution  (loading) 
of  each  blade  and  the  corresponding  pressure  field  generated  by  the 
mutual  Interference  of  the  blades  will  thus  be  obtained. 

In  Reference  5  Kemp  and  Sears  have  shown  that  the  velocity# 

V  #  Induced  at  a  rotor  blade  by  the  motion  of  the  steady  stator 
circulation  is  given  by 
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for  the  case  of  a  flat-plate 
distr  ibut  ion 

for  the  case  of  an  elliptic 
circulation  distribution 


a  ■ 


X 


s 


V  ■ 

r 


b  ■ 

P 


a  " 


U  - 

V  - 

r,s 


2c  2c 

X*  s 
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rotor  and  stator,  respectively. 

+  i(n/2-a  ) 

Tfae 

S 

circular  frequency  of  motion  ■  — ^ 

axial  distance  between  the  propeller  planes 

angle  of  blade  stagger  (the  complement  of  the  effective 
pitch  angle ) 

rotational  speed  of  rotor,  positive  down 
Inflow  speed  riative  to  blade 


Bessel  function  of  zero  and  first  order,  respectively. 
x,y*  Cartesian  coordinates  of  point  in  space  referred  to  a 

system  of  axes  fixed  in  the  blade  with  the  origin  at  the 
center  (see  Figure  2b)  and  with  positive  x  in  the 
direction  of  the  trailing  edge. 

As  mentioned  before,  the  additional  effects  on  the  rotor  blades  of  the 
vortex  wake  shed  by  the  stator  blades  are  also  taken  into  account  in 
the  above  reference.  In  this  case  the  Induced  velocity  v 
indicates  the  effect  of  wake  vorticlty)  at  the  rotor  blades  is  given 
in  complex  Fourier  series  form  byi 


t 

Z  Z  \k  tiv^.k(t-  x/v^)] 
m  k 


(18) 


where 

r 

r 

0 _ 

’'mk  d  cos  a 
s  s 


exp  [-Imrrc 


b  U 

'v  c  cos  a 
s  r  j 


r 


V  cl  . 
m 

r  s 


XT 


l+cot' 


r  s 


u 


V  c  /V 
r  r  r 


reduced  frequency 
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for  a  flat  plate  circulation  distribution 


for  an  elliptic  circulation  distribution 


^  =  summation  over  the  number  of  blades.  Theoretically  m  varies 

m  from  one  to  infinity,  however  only  a  few  terms  of  the  series 

^2 

contribute  to  the  induced  velocity  v 


L  =*  summation  over  k  from  -  ao  to  qd  omitting  k  =  0  , 
This  summation  actually  indicates  that  the  vortex  wake 
shed  by  the  stator  is  made  up  of  the  superposition  of 
waves  of  all  harmonics  from  -  qd  to  +  oo  . 

5 

Finally  the  expression  of  the  velocity  v  at  the  stator  blade  Induced 


by  the  motion  of  the  steady  rotor  circulation  F 

o 

Reference  5  as 


is  given  in 
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where  x^  is  referred  to  a  system  of  axes  fixed  in  the  stator  blade 
positive  toward  the  trailing  edge,  2TTU/d^  and  the  other  symbols 

are  as  defined  previously. 

The  foregoing  expressions  for  the  induced  velocities  at  the 
rotor  blades  and  stator  blades  can  be  written  in  a  more  convenient 
form  by  replacing  the  exponentials  with  well-known  series.  For 
instance,  the  exponent  which  appears  in  eq,  (17)  can  be  written  as 
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With  the  transformation 
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the  corresponding  exponential  is  written  as 
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Therefore,  the  unsteady  induced  velocity  on  a  rotor  blade  due  to  the  steady 
state  circulation,  i.e.  eq,  (17)  restricted  to  the  time-dependent  term 
alone^will  be  given  by 
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where  all  symbols  are  as  defined  in  eq.  (17)  and  eq,  (20). 

In  a  similar  fashion  the  unsteady  Induced  velocity  at  a  rotor 
blade  due  to  the  vortex  wake  shed  by  the  stator  blades  (see  eq,  18)  and 
the  velocity  on  the  stator  blade  induced  by  the  motion  of  the  steady 
rotor  circulation  will  be  given  by 


iv  kt 


'  ■  I  Z  ’mk'  '  ^  Z 


m  k 


and 


-  20 


s  iv  mt  m 

y  G  e  ® 

/ 1  in  os 

m 


OD 


z 


n  m 

(-1)  J^(i.Lg)cos  n0] 


(23) 


respectively,  where 
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In  References  6  and  7  the  acceleration  potential  of  a  thin  airfoil 
oscillating  with  an  arbitrary  mode  or  with  an  arbitrary  velocity  distri¬ 
bution  has  been  derived  in  the  X.  -plane,  in  which  plane  a  line  segment 
of  the  physical  z-plane,  from  x  •»  -1  to  x  »  1  ,  has  been  transformed 
into  a  circle  of  unit  radius  as  shown  in  Sketch  1  • 

For  an  airfoil  with  an  arbitrary  velocity  distribution 
io)t  n 

v(x,0)  -  W  e  (P^+  2  ^(-1)  )  (25) 


the  corresponding  acceleration  potential  is  given  by 


2UW[C(k)(P  +P,  )-P- ]  - ^  +  2UW  y  (4^  P  +p  -  p  )£i2j£  (26) 
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where  C(k)  •  Theodorsen's  function  of  argument  k 


K^(ik) 


K, (ik)  and  K  (ik)  ■  Modified  Bessel  functions  of  the  second 
i  o 

kind  of  order  1  and  zero,  respectively 
k  ■  reduced  frequency 

U  ■  velocity  of  incoming  flow  in  the  direction  of  the  airfoil  chord. 


■  polar  coordinates  of  the  field  point  from  the  leading 
edge,  (see  Sketch  1) 

r,0  ■  polar  coordinates  of  the  field  point  from  the  center  of 

the  airfoil. (see  Sketch  l) 


(it  must  be  noted  that  all  linear  dimensions  in  the  above  formulas 
are  expressed  in  terms  of  the  semi-chord  c  of  the  foil.) 


By  the  definition  of  the  acceleration 
potential  the  pressure  signal  at  any 
point  in  the  ^-plane  will  be  given  by 

P  -  -  pSi(  (27) 

and  the  lift  distribution  on  the  foil 
will  be  obtained  from  eq.  (27)  by  putting 
r-1,  0^  in  eq,  (26).  Now  if 

eq.  (25)  is  compared  with  the  equations 
for  the  Induced  velocities  on  rotor  or 
stator,  it  is  easy  to  identify  the 
corresponding  terms  in  each  case.  For 


SKETCH  I 


the  case  of  flow  (r^^)  of  the  steady  stator  circulation  past  the  rotor, 
the  various  terms  can  be  Identified  as  follows i 
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It  is  easy  to  show  that 


■  0  for  all  n  ^  1 

Hence  the  corresponding  acceleration  potential  in  the  ^i-plane  will  be 
given  by 
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The  unsteady  pressure  signal  in  non-dimensional  form  P  /j, 

1^ 

at  a  point  with  polar  coordinates  r^^,  0^^  in  the  ^-plane  will  be 
given  by 
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where  the  propeller  disc  loading  coefficient 
^  mpS2 
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has  been  used.  In  a  similar  fashion  the  pressure  signals  due  to  flow 
(r2)i  which  concerns  the  effects  of  the  stator  wakes  on  the  rotor,  and 
to  flow  (s)  of  steady  rotor  circulation  on  the  stator  will  be  given 
by  the  following  forms,  respectively 
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where  all  the  linear  dimensions  outside  the  first  parenthesis  are  expressed 

as  fractions  of  the  semi-chord  at  0.7b.  It  must  be  kept  in  mind  that 

P  and  P  are  referred  to  coordinate  systems  fixed  at  rotor  and  stator 
r  s 

foils  respectively.  The  function  sine^/r^^  is  determined  for  points 
of  known  polar  coordinates  in  the  ^-plane  of  each  foil.  It  can  be 
generally  stated  that  the  function  sine^/r^  is  the  field  point  function 
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which  depends  also  on  the  geometric  characteristics  of  the  propeller 
and  the  remaining  terms  constitute  the  loading  function. 

For  computational  simplification  it  will  be  assumed  that  at  time 
t  =  0  the  center  of  the  chord  of  stator  or  rotor  lies  on  the  x-axis 
of  the  coordinate  system  fixed  relative  to  ship  and  the  counter¬ 
rotating  propeller.  Then  by  a  simple  trigonometric  transformation  a 

point  with  coordinates  x  ,  y  in  the  original  system  will  be  determined 
^  0  0 

with  respect  to  the  orthogonal  axe's  fixed  in  the  rotor  or  stator  blades 
as  followsi 

For  instance  for  the  case  when  the  origin  of  the  blade  axes  is  on  the 
x-axis  of  the  ship  the  coordinates  of  a  point  will  be 

X  =  (x  -b  /2)cosa  +  y  sin  a 

y  =  -(x  -b  /2)sina  +  y  cos  a 
op  T  O  T 

with  respect  to  a  system  fixed  in  the  rotor,  and 

X  ■  (x  +b  /2)cos  a  -  y  sin  a 
sop  s  o  s 

y  =  (x  +b  /2)sin  a  -  y  cos  a 

‘^S  Op^  3*^0  S 

with  respect  to  the  stator,  where  b^  is  the  distance  between  the  two 

propeller  planes.  Since  the  pressure  signal  is  determined  by  means  of 

eq,  28,  29  and  30  in  terms  of  sin0j^/r^  (i.e.  in  the  transform  ^-plane), 

the  transformation  from  the  physical  plane  (x  ,y  )  or  (x  ,y  )  is  the 

X*  X*  S  5 

next  essential  step.  For  a  more  systematic  evaluation  of  the  field 

point  function  sin9,/r,  its  contours  for  various  values  of  (x  ,  y  ) 

^  1'  1  .i^s  “^rs 

have  been  plotted  in  Fig.  (3).  This  chart  is  derived  from  the  Joukowskl 
transformation 

K  ^  z  1 

where  z  =  x+iy,  K  ^  S+lq  ,  and  x  and  y  are  expressed  in  terms 
of  semi-chord  at  the  effective  radius.  The  following  relations  have 
been  used  in  evaluating  the  contours  of  sin9^/r^ 


FT 


5  +71 


Ti  -  sin  0^ 


K  =»  r^cos0^  -  1 

The  pressure  signal  due  to  interference  effects  between  the  two  propellers 
can  now  be  determined  by  means  of  eq.  28,  29  and  30  and  by  the  contours 
of  the  field  function  sin0^/r^  . 

Numerical  calculations  are  made  of  the  pressures  along  four 
longitudinal  transverses  at  various  axial  distances  for  the  propeller 
configuration  3-0-3  (no.  2714,  see  ref,  8)  with  the  following  pertinent 
characteristics 


3-0-3  (No.  2714) 


No.  of  blades,  m  =  3 
Diameter,  2b  •>  0,68  ft. 
Semi-chord,  c  ,  at  0.7b 
Distance  between 

propeller  planes,  b^ 
RPM  =  900 


0.186  ft. 

0.124  ft. 


Speed  of  advance  =  8,16  ft/sec. 


Results  of  these  calculations  are  presented  in  Tables  I  and  II. 

Table  I  gives  the  blade  loading  factor  (ratio  of  unsteady  dimensionless 

pressure  to  the  field  point  function)  for  rotor  and  stator  blades. 

Table  II  gives  the  field  point  function  referred  to  rotor  and  stator 

blades  for  one  position  of  the  effective  blade  chord.  From  the  products 

of  stator  or  rotor  blade-loading  and  field-point  factors  for  a  given 

x^,  y^  are  obtained  the  steady  and  unsteady  dimensionless  pressure  signals 

at  that  point,  per  blade  of  stator  or  rotor.  Figure  4  shows  the  sum 

of  the  blade-frequency  pressure  signals  due  to  rotor  and  stator  blades 

plotted  versus  axial  distance  different  tip  clearances 

r/b“  unsteady  interference  effects  two  circulation 

distributions  are  assumed,  flat  plate  and  elliptic.  The  results  of  Fig,  4 

have  been  obtained  for  the  blade  chord  center  arbitrarily  placed  on  the 

X  -axis  (y^-  O), 

0  0 


-  25 


If  the  pressure  signal  due  to  the  interference  effects  is  compared  with 
the  quasi-steady  pressures  computed  by  the  method  of  Part  I  (see  Table 
III)  for  the  counterrotating  system  treated  as  though  the  two  propellers 
operate  as  separate  units  in  three-dimensional  flow,  it  is  seen  that 
up  to  a  distance  of  0.4  propeller  radius,  where  the  quasi-steady  pressure 
(three-dimensional)  attains  its  maximum,  that  pressure  is  at  least 
10  times  the  pressure  signal  due  to  the  interference  contribution.  At 
larger  distances,  however,  no  conclusion  can  be  drawn  as  to  the  relative 
magnitudes  of  the  pressure  signal  contributed  by  the  interference  effect 
and  of  the  quasi-steady  signal  because  the  rate  of  attenuation  with 
axial  distance  of  the  two-dimensional  pressure  is  considerably  smaller 
than  the  rate  of  decay  of  the  three-dimensional  one.  It  can  be  stated, 
however,  that  in  the  actual  three-dimensional  case  the  interference 
effects  will  be  smaller  and  will  decay  more  rapidly. 

To  decide  more  conclusively  about  the  relative  importance  of 
the  interference  effect  a  two-dimensional  analogue  of  the  uncoupled 
propeller  configuration  is  considered.  Expressions  for  the  blade 
frequency  pressure  signal  contributed  by  the  forward  and  after  string 
of  blades  operating  independently  have  been  developed  by  W,  Jacobs 
by  a  procedure  similar  to  that  used  in  Reference  5  to  determine  the 
corresponding  pressure  due  to  the  interference  effects. 


In  Appendix  D  it  is  shown  that  the  dimensionless  pressure 
signals  at  a  field  point  due  to  the  motion  of  cascades  of  blades  of 
the  fqrward  and  after  propeller  operating  independently  are  given  as 
follows  I 


(1) 


4TtVj^  sin9. 


nS.  io  ^imvt  r 

^  Re  ^  e  e  |c(p)[j^(p)-lJj(p)li-U^(p)| 


(31) 


ahd 

(2) 


4nV  sin©  ia  +lmvt  r 

X  ®  ®  [lp(X)“Il(X)]  jc(p)[J^(p)-iJ^(p)]*iJ^(p)j. 

(32) 


respectively, 

where 

nig  ■  number  of  blades 

SJ  “  angular  velocity  “  2rt  (r.p.s.) 

V_  ■  resultant  velocity  along  the  blade  chord 
R 

i  o 

-  2nmce 

- a— 

and  all  geometric  characteristics  are  measured  at  0.7  propeller  radius 
and  are  defined  previously.  Figure  5  graphically  exhibits  the  comparison 
in  a  two-dimensional  mathematical  model  of  the  total  blade  frequency 
pressure  signal  obtained  from  the  motion  of  the  two  rows  of  blades 
operating  separately  with  the  pressure  generated  from  the  interference 
effects.  It  is  seen  from  this  figure  that  the  propeller  interference 
gives  rise  to  a  pressure  signal  at  most  15%  of  the  corresponding  pressure 
obtained  from  the  summation  of  the  pressure  fields  of  the  two-rows  of 
cascade  blades  operating  as  separate  units  (uncoupled).  From  the  results 
it  can  be  assumed,  for  practical  purposes,  that  the  additional  pressure 
contributed  by  the  mutual  blade  interference  effects  can  be  neglected. 

The  limited  numerical  results  presented  here  indicate  thati 

1,  the  effect  of  the  vortex  wakes  shed  by  the  stator  blades  on 
the  rotor  (aft  propeller)  is  comparable  in  magnitude  to  the  effect 
induced  by  the  stator  (forward  propeller)  blades  themselves, 

2,  the  contribution  of  the  rotor  to  the  stator  is  small, 

3,  the  interference  part  of  the  vibratory  pressure  signal  is 
a  small  part  of  the  total  pressure  signal  in  the  near  hydrodynamic 
field  and  can  be  neglected  for  the  purpose  of  this  analysis. 
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Part  III 

Calculation  of  the  Pressure  Field 
Near  3-0-3  Counterrotating  Propellers 

The  vibratory  pressure  of  various  points  in  the  neighborhood  of 
a  3-0-3  counterrotating  propeller  is  determined  by  utilizing  the  two 
expressions  evolved  in  Part  I  for  the  pressure  signal  when  the  system 
of  dual  propellers  is  "in-phase"  and  "out-of-phase". 

In  figure  (2a)  the  two  3-bladed  propellers  are  shown  to  be 
located  at  with  respect  to  the  origin  of  the  coordinate  system 

fixed  on  their  common  axis  of  rotation  at  the  point  midway  between  the 
propeller  planes.  In  the  computational  work  the  two  counterrotating 
components  are  No.  2714  propellers  (Reference  8)  whose  pertinent 
characteristics  have  been  given  previously  in  Part  II,  The  distance 
between  them  is  taken  as  0,364  of  the  propeller  radius,  l,e,  ^0.182 
in  radii. 

The  locations  at  which  the  pressure  signal  is  'jvaluated  are  given 
in  the  same  figure.  The  pressure  is  computed  along  four  longitudinal 
transverses  having  tip  clearances  t/d  »  0,1,  0,15,  0,20,  0,30  and  at 
6  different  stations  along  each  layer,  having  axial  distancesi 
0.182,  0,382,  0.482,  0.582,  and  0.782  of  the  propeller  radius. 

Expressions  11,  12,  14  and  15  which  have  been  developed  in  Part  I 
can  be  used  in  evaluating  the  harmonic  coefficient  of  each  individual 
propeller  with  special  precaution  being  taken  for  the  axial  distance  x  , 
For  the  forward  propeller  the  x  is  replaced  by  snd  for  the 

after  propeller  by  x+x^  ,  since  the  general  expressions  have  been 
derived  for  the  origin  of  the  coordinate  axes  at  a  propeller  center. 

For  several  stations  the  coefficients  of  the  first  and  second 
blade-frequencies  were  evaluated  by  the  developed  expressions  and  their 
values  were  found  to  be  very  close  to  those  computed  by  numerical 
methods.  However,  as  was  mentioned  before,  the  computational  work  was 
donein  entirety  on  the  ELECOM-100  computer,  since  the  coding  was  available. 
The  total  pressure  signal  was  determined  by  summing  up  the  contributions 
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of  three  individual  blades  located  at  120°  intervals.  Then  the  harmonic 
analysis  was  performed  by  the  convenient  numerical  procedure  described 
in  Reference  4, 

After  the  coefficients  of  the  first  and  second  harmonics  are 
determined,  the  pressuru  signal  in  the  field  of  the  3-0-3  counterrotating 
propeller  is  evaluated  by  using  eq,  (9)  and  (lO),  when  blades  are 
"in-phase"  and  "out-of-phase",  respectively. 

Results  for  the  total  anplitude  of  the  first  and  second  harmonics 
for  the  above  two  conditions  have  been  plotted  in  Figures  6  and  7  versus 
tip  clearances,  at  various  axial  distances. 


Cross-plots  of  these  results  are  shown  in  Figures  8,  9  and  10, 
where  in  addition  the  corresponding  amplitudes  of  a  single  3-bladed 
propeller  producing  equivalent  thrust  and  having  the  same  RPM  have  been 
plotted  for  comparison.  The  diameter  of  the  conventional  propeller  has 
been  determined  by  calculations  made  at  theDavld  Taylor  Model  Basin  on 
the  basis  of  the  Troost  optimum  criteria,  for  SHP  “  15,000  at  200  RPM, 

It  is  found  on  this  basis  that  the  diameter  of  an  equivalent 
single  propeller  is  given  in  terms  of  the  propeller  diameter  of  the 
counterrotating  system  as 

°s  =  ^CR 

Therefore  the  corresponding  thrust  loading  coefficient  T'  of  the  single 

6 

propeller  will  be  given  by 


10^(125)2  ^^-25)' 


T»  •  1.28  T' 
XR  ‘CR 


The  amplitudes  of  the  pressure  signal  generated  by  the  equivalent  con¬ 
ventional  propeller  are  evaluated  at  the  same  field  points  at  which 
the  pressure  signals  of  the  counterrotating  system  have  been  determined. 
The  results  of  these  calculations  are  presented  in  Tables  III  and  IV 
as  well  as  graphically  in  Figurffi  8,  9  and  10,  This  comparison  indicates 
that  for  the  first  harmonic  when  the  blades  are  in  alignment  the  vibratory 
characteristics  of  the  counterrotating  configuration  are  undoubtedly 
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superior  to  the  corresponding  characteristics  of  the  equivalent  con¬ 
ventional  single  propeller  in  the  immediate  neighborhood  of  the  propeller 
i,e,  x/b  <  0.6  and  t/d  <  0,3.  Outside  of  that  region,  however,  (see 
Fig,  8  and  Table  III)  the  equivalent  single  propeller  is  more  advantageous 
since  the  counterrotating  configuration  in  the  "in-phase"  condition  acts 
as  a  single  propeller  with  the  same  number  of  blades  as  one  stage  of  the 
dual.  In  the  out-of-phase  situation,  however,  the  counterrotating 
propeller  system  has  always  better  vibratory  characteristics  than  the 
equivalent  single  propeller  as  shown  in  Figures  9  and  10  and  Tables  III 
and  IV,  Since  the  in-phase  situation  is  relatively  infrequent,  it  can 
be  generally  stated  that  the  counterrotating  propeller  configuration 
is  superior  to  an  equivalent  conventional  single  propeller  as  far  as 
the  vibratory  characteristics  are  concerned. 

Of  course,  mechanical  complications  are  involved  in  the  counter¬ 
rotating  propeller  configuration  which  are  not  overlooked  in  selecting 
the  most  advantageous  configuration.  It  must  also  be  pointed  out,  as 
was  mentioned  in  the  Introduction,  that  other  criteria,  like  the  avoidance 
of  cavitation,  noise  reduction,  increase  in  efficiency,  should  be  con¬ 
sidered  together  with  the  vibratory  characteristics  in  judging  whether 
or  not  the  advantages  that  are  indicated  for  a  counterrotating  system 
counterbalance  the  disadvantages. 

It  is  interesting  to  notice  from  the  polar  diagram  (Figure  11)  of 
the  first  blade  harmonics  of  the  pressure  signal  as  a  function  of  the 
observer's  position  relative  to  blade  overlap,  that  the  pressure 
at  the  lines  of  blade  intersections  (in-phase)  has  the  frequency  and 
the  mode  of  a  single  propeller  having  the  same  number  of  blades  as 
one  stage  of  the  dual,  whereas  in  the  direction  midway  between  the 
Intersection  lines  (out-of-phase)  the  pressure  is  that  of  a  single 
propeller  having  the  total  number  of  blades  of  the  set  of  two.  In 
this  figure  the  polar  diagram  of  the  pressure  is  made  up  of  six  loops 
since  there  are  six  overlaps  of  the  counterrotating  blades  in  one 
revolution.  In  the  out-of-phase  condition,  there  are  six  directions 
which  coincide  with  the  bisectors  of  the  60°  intervals  between  overlaps. 


CONCLUSIONS 


This  investigation  indicates  that  the  counterrotating  propeller 
system  has  vibratory  characteristics  much  superior  to  the  equivalent 
single  propeller  having  the  same  number  of  blades  as  one  of  its 
components^  delivering  the  same  thrust  at  the  same  RPM  at  optimum 
conditions.  At  a  distance  approximately  0.4  -  0,5  of  the  propeller 
diameter  from  the  propeller  plane  the  conventional  propeller  becomes 
better  than  the  in-phase  counterrotating  system,  but  at  Instances  when 
the  blades  are  not  in  coincidence,  a  situation  which  occurs  more  often 
in  a  complete  revolution,  the  superiority  of  the  3-03  system  Is 
evident  at  points  farther  away. 

This  study  indicates  also  that  the  vibratory  pressure  signal  of 
the  counterrotating  system  can  be  obtained  from  the  summation  of  the 
pressure  fields  of  its  two  components,  as  though  the  two  propellers 
operate  as  separate  units,  since  their  mutual  interference  effects 
calculated  on  the  basis  of  2-dimensional  theory  have  been  found  to  be 
a  small  part  of  the  total  pressure  signal.  In  an  actual  three- 
dimensional  situation,  the  contribution  of  the  interference  effect 
to  the  vibratory  pressure  signal  will  be  smaller. 

The  study  of  the  blade-blade  interference  effects  in  the  two- 
dimensional  analogue  indicates  that  the  effect  on  the  aft  propeller  of 
the  forward  propeller  with  its  shed  vorticities  is  much  more  Important 
than  the  effect  of  the  aft  propeller  on  the  forward  one,  hence  the 
latter  effect  can  be  omitted  as  negligible. 
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TABLE  I 


P  sinQj^ 

— )  per  blade  of  rotor  for  J  »  0.8 


Blade  Harmonics 

1 
2 
3 

1 
2 

P  sine, 

_  /  ( - )  per  blade  of  stator  for  J  =  0,8 

^1  s 

1)  Effect  of  Steady  circulation  of  rotor 

-.021  -  1.001 


With  Flat  Plate  Distribution  With  Ell IpticDistr Ibutlon 
_ of  Circulation _  of  Circulation _ 

1)  Effect  of  steady  circulation  of  stator 

-.057  +  i.005  -.049  +  1.138 

i.006  .002  +  i.019 

i.OOl  .001  +  1.006 

2)  Effect  of  Wake  Vorticities  Shed  by  Stator 

+.010  +  i.027  .022  +  1.010 

-.004  -  1.037  -.022  -  1.002 


1 

2 

3 


-.001  -  i.026 
-.001  -  i.008 
-.001  -  i.003 


002 


TABLE  II 


Field  point 

sinQj^ 

function  ( - 

’^l 

-)  referred  to  rotor  and 

stator  blades 

V* 

(«  r/R) 

(center  of 

X  /r 
o' 

blade  chord  on  x  -axis) 
0 

sinG 

1  r 

SinG, 

1  s 

1.2 

0 

.030 

.200 

.182 

.050 

.167 

.382 

.073 

.120 

.488 

.077 

.083 

.582 

.081 

,050 

.782 

.093 

-.005 

1.6 

.093 

-.080 

2.0 

.088 

-.080 

3.0 

.083 

-.070 

4.0 

.060 

-.055 

6.0 

.042 

-.039 

8.0 

.033 

-.030 

1.4 

0 

.023 

.150 

.182 

.038 

.133 

.382 

.055 

.100 

.488 

.062 

.073 

.582 

.067 

.050 

.782 

.080 

.013 

1.6 

.090 

-.060 

1.6 

0 

.019 

.130 

.182 

.033 

.111 

.382 

.044 

.080 

.488 

.050 

.067 

.582 

.055 

.050 

.782 

.067 

.020 

1.6 

.087 

-.040 

TABLE  III 


First  Blade  Harmonics  of  Pressure  Signal  pAqj^  due  to 

1)  3-0-3  counterrotating  propellers  of  radius  b 

2)  3-bladed  single  propeller  of  radius  1.25b 
(equivalent  thrust) 


Field  point 

posit  ion 

3-0-3  in-phase 

3-0-3  out-of-phase 

Single 

< 

O 

X 

r/b 

’’AiR. 

0 

1.2 

.083 

.135 

-.070 

1.3 

.060 

.072 

.071 

1.4 

.044 

.041 

.051 

1.6 

.026 

.015 

.029 

.182 

1.2 

.097 

.063 

.47 

1.3 

.068 

.041 

.23 

1.4 

.048 

.026 

.094 

1.6 

.027 

.012 

.035 

.382 

1.2 

.12 

.041 

.16 

1.3 

.078 

.019 

.12 

1.4 

.054 

.010 

.087 

1.6 

.028 

.006 

.041 

.488 

1.2 

.10 

.044 

.11 

1.3 

.071 

.024 

.092 

1.4 

.050 

.013 

.071 

1.6 

.027 

.005 

.036 

,582 

1.2 

.082 

.039 

.086 

1.3 

.061 

.024 

.072 

1.4 

.045 

.015 

.055 

1.6 

.026 

.006 

.032 

.782 

1.2 

.047 

.024 

.041 

1.3 

.039 

.017 

.035 

1.4 

.031 

.012 

.030 

1,6 

.020 

.006 

.020 

♦ 

Subscript  CR  refers  to  counterrotating  propellers 


TABLE  IV 


Second  Blade  Harmonics  of  Pressure  Signal  P/^T*  due  to 

UK 

1)  3-0-3  counterrotating  propeller  of  radius  b 

2)  3-bladed  single  propeller  of  radius  1,25b 
(equivalent  thrust) 


Field  point 

pos it  ion 

S-0-3  (in-phase  or  out-of-phase) 

•Single 

X  /b 

0 

r/b 

0 

1.2 

.032 

.027 

1.3 

.019 

.043 

1.4 

.011 

.028 

1.6 

.004 

.011 

.182 

1.2 

.035 

.23 

1.3 

.020 

.12 

1.4 

.011 

.056 

1.6 

.005 

.013 

.382 

1.2 

.041 

.054 

1.3 

.021 

.045 

1.4 

.012 

.032 

1.6 

.004 

.012 

.488 

1.2 

.028 

.035 

1.3 

.017 

.029 

1.4 

.010 

.021 

1.6 

.004 

.009 

.582 

■  1.2 

.018 

.021 

1.3 

.012 

.017 

1.4 

.008 

.012 

1.6 

.003 

.006 

.782 

1.2 

.007 

.007 

1.3 

.005 

.005 

1.4 

.004 

.004 

1.6 

.002 

.003 

♦ 

Subscr Ipt 

CR  refers  to 

counterrotating  propellers 
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FIGURE  4.  AMPLITUDE  OF  THE  lit  BLADE  HARMONIC  OF  THE  UNSTEADY  DIMENSIONLESS  PRESSURE  SIGNAL  P/T* 
PER  BLADE  DUE  TO  INTERFERENCE  EFFECTS,  FOR  A  TWO-DIMENSIONAL  MATHEMATICAL  MODEL 
REPRESENTING  THE  COUNTERROTATING  SYSTEM.  (Jt  0.8,  CENTER  OF  BLADE  CHORD  ON  X,- AXIS) 
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FIGURE  6.  AMPLITUDES  OF  THE  Isl  BLADE  HARMONIC  OF  THE  DIMENSIONLESS 
PRESSURE  IP/T'I  due  TO  THE  UNCOUPLED  COUNTERROTATING 
SYSTEM  FOR  THE  IN  AND  OUT- OF- PHASE  CONDITIONS 
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FIGURE  8.  COMPARISON  OF  AMPLITUDES  OF  THE  lit  BLADE  HARMONIC  OF  THE  DIMENSIONLESS 
PRESSURE  DUE  TO  THE  COUNTERROTATING  SYSTEM  IN  THE  IN-PHASE  CONDITION 
WITH  THOSE  DUE  TO  THE  EQUIVALENT  CONVENTIONAL  PROPELLER 


FIOURE  9.  COMPARISON  OP  AMPLITUDES  OF  THE  lit  BLADE  HARMONIC  OF  THE  DIMENSIONLESS 
PRESSURE  DUE  TO  THE  COUNTERROTATING  SYSTEM  IN  THE  OUT- OF- PHASE  CONDITION 
WITH  THOSE  DUE  TO  THE  EQUIVALENT  CONVENTIONAL  PROPELLER 
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FIGURE  10.  COMPARISON  OF  AMPLITUDES  OF  THE  2nd  BLADE  HARMONIC  OF  THE  DIMENSIONLESS 
PRESSURE  DUE  TO  THE  COUNTERROTATING  SYSTEM,  IN  AND  OUT-Of-PHASE,  WITH 
those  due  to  THE  EQUIVALENT  CONVENTIONAL  PROPELLER 


I80» 


165* 


I 

I 

I 

1 

I 


I98* 


FIGURE  M.  POLAR  DIAGRAM  OF  THE  Ift  BLADE  HARMONIC  OF  THE  DIMENSIONLESS  PRESSURE  DUE 
TO  THE  COUNTERROTATING  SYSTEM,  IN  THE  IN  AND  OUT-OF-PHASE  CONDITIONS,  AT 
r/b  •  1.2,  Xg/b  «  0.382 
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APPENDIX  A 


It  is  known  that  the  sum  of  a  finite  geometric  series, 

o  1  n+1 

^  2^  n  1-z 
1+z+z  +...,z  “ 


1-z 


iS 

If  z  ■  e  ,  it  can  be  shown  that 
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These  identities  are  used  to  prove  that  the  vibratory  pressure  is  made 
up  of  all  the  harmonics  whose  order  is  an  integer  multiple  of  the  number 
of  blades. 

The  cosine  part  of  the  pressure  distribution  is  given  by 
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where  the  first  term  is  due  to  the  first  blade  and  the  second  is  due  to  the 
remaining  (m-l)  blades. 
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Proof  1 
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APPENDIX  B 


Expressions  for  the  pressure  signal  due  to  a  system  of  counterrotating 

propellers  will  be  developed  at  the  instances  when  the  blades  are  in 

coincidence,  "in  phase",  and  In  the  "out-of-phase"  position.  For  the 

"in-phase"  condition  t  ■  —  ,  whereas  the  out-of-phase  situation  occurs 
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APPENDIX  C 


The  derivation  of  the  Fourier  coefficients  of  the  second  blade 
frequency  of  the  pressure  signal  near  a  three-bladed  propeller  given  by 
eq.  (14)  and  (15)  will  be  shown  in  very  general  steps  and  the  most 
representative  type  of  integration  will  be  carried  out. 

After  the  designated  integration  with  respect  to  s  the  instantaneous 
pressure  signal  given  by  eqs.  (l:)  and  (2)  becomes 
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To  evaluate  the  second  harmonic  of  expression  (C-l)  it  is  necessary  to 
find  the  sine  and  cosine  Fourier  coefficients.  The  cosine  coefficients  arei 
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coefficients  of  the  components  (l^^,  I^, . I^)  of  eq,  (C-l) 

On  expanding  cosda  and  sindo  in  terms  of  cosa  and  sina  the 
following  representative  types  of  integral  will  occuri 
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where  all  A's  can  be  determined  by  the  recurrence  fornula  of  eq.  (14) 
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The  E  integral  will  be  shown  to  be  a  sum  of  all  three  kinds  of  complete 
in 

elliptic  integrals.  The  technique  of  transforming  this  Integral  into  the 
complete  elliptic  integrals  is  given  in  Reference  2, 
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In  Reference  10,  by  the  same  general  procedure,  the  various  I's  were 
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where  are  given  in  eq,  (14)  and  A's  by  the  recurrence  formula 

(see  eq,  14),  The  results  have  been  written  in  rathar  extended  form,  in 
order  to  show  the  structure  of  the  different  expressions  so  that  general 
rules  may  be  set  down  for  the  recapture  of  an  expression  from  the 
preceding  one. 

The  following  generalities  can  be  deduced  for  the  above  express  Ions i 
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5,  The  highest  subscript  of  A  in  a  given  expansion  Is  four 

less  than  the  exponent  of  the  original  integral  i,e,  (n-4)  so 

that  the  A  terms  are  from  A  to  A  .  , 
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APPENDIX  D 
by  W,  Jacobs 

This  appendix  describes  the  use  of  the  two-dimensional  mathematical 
model  of  the  counterrotating  propeller  configuration  in  order  to  obtain 
the  pressure  signals  due  to  the  independent  motion  of  each  of  the  two 
cascades  of  blades  representing  the  propellers. 

The  coordinate  axes  are  as  in  Sketch  2, 


If  the  field  point  (x,y)  is  referred  to  the  axes  in  a  blade  of  an  infinite 
cascade*  then  its  velocity  potential*  referred  to  the  singularities  (C*0) 
on  the  blade*  is  given  by 
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where  y  =  the  strength  of  the  circulation 
Ti  -  Ut 

U  “  rotational  speed  of  the  propeller,  ft/sec. 
n  =  the  blade  number  in  the  infinite  cascade  and  the  other 
symbols  are  as  defined  before. 
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There  are  two  possibilities: 
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The  velocity  distribution  along  the  blade  chord  will  be  obtained  by  setting 
y*0  ,  X  “  c  COS0  •  Then 
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By  using  the  acceleration  potential,  the  pressure  can  be  obtained 
as  before,  at  any  point  defined  by  the  field  point  function  (sinSj^/r^  • 
The  dimensionless  pressure  signal  due  to  the  velocity  distribution  (6) 
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where  m^  ■  number  of  blades  in  the  propeller 
a  -  2n (r .p.s. ) 

V  -  resultant  velocity  along  the  blade  chord 
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and  all  geometric  characteristics  are  measured  at  0.7  propeller  radius. 

2)  (5-x)coso  +  y  sino  >  0  in  which  case  the  time-dependent  part  of  eq.  (4) 


becomes 


V  -  Re 

ncd 


!{•“ 


2Ttme^°x  -12"me^%  .  _2Trme^5 

e  ^  e  ^  e~  ^  f /c^  e 


D-5 


Let  5  ■  -c  cosQ  ,  then  dK  ■  csinQdG  and  when  5  •  -c  ,  0“O  and  5  ■  , 

0“«  ,  Then 


-2Ttnie^°? 

d? 


Tt 

c  /  (l+cos6)e 
o 


2ttme^°c  COS0 


d0 


cn 


{T  /2Trmce^^  ^  _  /2nmce^v\ 
o''  d  *  1  ^  d 


The  velocity  distribution  along  the  blade  chord  will  be  obtained 
by  setting  y  *  0  ,  x  =  -c  cos0  .  Then 

r  ^  f  la  +  imvt  ^  la  ') 

v(x,0>-Sg  Re  ^  |e  e  exp  (-  )[ I^(X)*I^(X)]  | 

m 


or 


v(x,0)* 


Imvt 


0) 


m 


[l^(X)+I^(X)][j^(p)+2  Jjj(p)co8n0l 

n"l 


} 


(D-8) 


I 


and 


P 

V"  “ 


4TtVD  sin0 
in^  S2  d'  r 


lI^)Re^|e^e“^'’^[l^(X)*I^(X)l  {c(p)[j^(p)-lJj(n)]  ♦IJ^(p)}  | 

^  m 


(D-9) 


Whether  equation  (D-V)  or  equation  (D-9)  Is  used  to  calculater  the 
pressure  signal  at  a  point  x«y  referred  to  particular  blade-chord  axes 
depends  on  the  condition 

l)  (x-c  cosO)  cosa  ^  y  slno>  0 
or 


I 


2)  (-C  COS0  -  x)  cosii  +  y  sina  >  0 

at  that  point.  It  is  found  that,  for  the  case  of  blade  chord  centered 
on  the  x^-axis  and  for  points  on  the  ship's  center, plane  at  r/h  ■1.2 
and  x^/b  from  0  to  8  ,  eq.  (D-7)  applies  to  the  after  propeller  and 
eq,  to  the  forward  propeller. 
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